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We have studied the thickness-induced superconductor-to- insulator transition in the presence of a 
magnetic field for a-NbSi thin films. Analyzing the critical behavior of this system within the "dirty 
boson model", we have found a critical exponents product of ^ 0.4. The corresponding phase 
diagram in the {H,d) plane is inferred. This small exponent product as well as the non-universal 
value of the critical resistance found at the transition call for further investigations in order to 
thoroughly understand these transitions. 
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I. INTRODUCTION 

Low temperature transport in disordered conducting 
materials imply quantum interferences, Coulomb repul- 
sion, and superconducting fluctuations. Since 2D is the 
lower critical dimension for the existence of both the su- 
perconducting and the metallic states, transport proper- 
ties of such disordered thin films have attracted continu- 
ous attention since the 1960s in order to understand what 
ground states are allowed in those systems and study the 
nature of the quantum phase transitions between the dif- 
ferent phases 0, @, Q . 

Quantum Phase Transitions (QPT) occur when a pa- 
rameter in the Hamiltonian is varied, resulting in a 
change of the system's ground state. These transitions 
therefore take place at zero temperature and are driven 
by quantum fluctuations, contrary to classical phase 
transitions which are controlled by thermal fluctuations. 
Near a QPT, the quantum fluctuations have a charac- 
teristic lengthscale - the correlation length ^ - diverging 
as ^ (X where u is the correlation length critical ex- 
ponent, Sk = ^^^^""^ the distance of the considered sys- 
tem to the K-diWen transition, and K an experimentally 
tunable parameter which critical value is Kc. The fluc- 
tuations are also characterized by a vanishing frequency 
Q (X where z is the dynamical critical exponent. The 
two critical exponents u and z define the universality 
class to which the transition belongs. 

In the case of Superconductor to Insulator Transitions 
(SITs) in disordered thin films, the tunable parameter in 
the Hamiltonian can be the disorder or the magnetic field 
H. The most popular theoretical model to explain these 
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SITs is the "dirty boson model" developed by M.P.A. 
Fisher [2]. In this model, the coherence of the super- 
conducting state is destroyed by quantum fluctuations 
of the order parameter's phase and the system amounts 
to interacting bosons in the presence of disorder. The 
superconducting and insulating phases are then dual to 
one another : the superconducting phase consists of lo- 
calized vortices and condensed Cooper pairs, whereas the 
insulating phase is characterized by condensed vortices 
and localized Cooper pairs. Both disorder and magnetic 
field-driven transitions have similar description within 
this frame : in the quantum regime, for DC measure- 
ments, the sheet resistance obeys a scaling law that is 
solely dependent on the variable S * T~ ^ [jj, l2| : 

i?(5,T) = i?e/MT-^) (1) 

where is the critical sheet resistance and / an uni- 
versal scaling function having an unique constraint : 
/(O) = 1. a is a non universal constant 2; = 1 is ex- 
pected due to the long-range Coulomb interactions and 
the "dirty boson model" predicts > | = 1 as well as 
an universal value of the system's sheet resistance at the 
transition Rc = Rq = = 6500 Q [3] . Despite obeying 
to the same scaling laws (equation ([1])), the field-induced 
transition and the disorder-induced transitions have dif- 
ferent physical grounds : in the magnetic field-induced 
SIT, the vortex density increases with the magnetic field, 
until they delocalize and Bose condense ; in the disorder- 
induced SIT at zero field, the Bose condensation is un- 
dergone by the vortex/ ant ivortex pairs. These two SITs 
hence have no reason to have the same critical exponents 

Experimentally, number of disordered superconducting 
films experience a SIT when submitted to a perpendic- 
ular magnetic field. However, they do not all behave in 
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the same way. Following Gantmakher's comment [s*], one 
can separate them into two different categories. Some 
compounds exhibit an insulating phase which low tem- 
perature resistance is only 10% above their high tem- 
perature resistance. This behavior resembles more the 
one of a conductor in the presence of weak localiza- 
tion than the one of an actual insulator [6]. This is 
the case of Mo^Gei_^ [7], Mo^Sii-^ [8], Be @, a-Bi 
or Nd2-.rCe.^Cu04+7y [5|. Other systems, such as 
amorphous indium oxyde or TiN [l2|, have, in the 
same conditions, a much more important increase in re- 
sistance - up to a factor 10. Their resistance then have 
an exponential increase with the temperature [ll|, fl3l |. 
The renormalization analysis of these field-induced SIT 
gives 0.75 < vhz < 1.35, independently of the above- 
mentioned categories. 

The experimental realizations of the thickness-induced 
SIT, where tuning the system's thickness is taken to be 
a mean of varying its disorder, are far more rare because 
of the experimental difficulty of synthesizing microscopi- 
cally identical films which only differ by their thicknesses. 
In the case of this transition, the distinction previously 
made does no longer exist : all studied compounds show 
a drastic increase in resistance of many orders in mag- 
nitude when their thickness is lowered [lO|. However, 
one can make another distinction. Some systems, such 
as a-Bi [10], are very sensitive to any thickness variation 
: a fraction of angstrom difference engenders resistance 
increases of several orders of magnitude at low tempera- 
ture. This behavior is comparable to the one observed in 
granular systems [14]. On the other hand, systems such 
as MoC present a more progressive thickness-dependence 
[15.] . Values of the critical exponents have only been re- 
ported for a-Bi [io| : ~ 1.3. 

Whichever the parameter tuned to induce the SIT, and 
contrary to the predictions of the "dirty boson model", 
experiments show an important variation in the values 
of the critical sheet resistance at the transition Rc 0, B 
B [lO, [HI, [HI • Within one system, Rc can vary between 
2000 Vt to 9000 Vt [9] depending on the applied magnetic 
field or the normal resistance of the sample. Theories 
introducing a fermionic channel of electronic conduction 
have been developed to explain the non-universality of 
Rc [7] but these are not entirely satisfactory since they 
do not account for values of Rc larger than Rq [l^. 

As one can see, all the experimental realizations of the 
SITs in thin disordered films show a large variation in 
the measured critical exponents, as well as in the critical 
resistance. This has led to the questioning of the "dirty 
boson model" . Some have suggested a percolation-based 
mechanism [l^, others the contribution of fermions to 
the conduction near the transition 0. Moreover, the 
fiat R{T) curves found near the transition have put into 
question Fisher's picture of an unique metallic separa- 
trix between the superconducting and insulating regimes. 
Some have suggested the existence of an intermediate 
metallic phase - the Bose metal. 

In this context, it seemed to us particularly interest- 



ing to provide another example of such transition. 2D 
Nba^Sii-a; films are interesting systems for this study. We 
have previously shown that these films experience a mag- 
netic field-tuned SIT [19] with a product of critical expo- 
nents vhz = 0.67, in agreement with other experimen- 
tal data [iO] but in contradiction with the "dirty boson 
model". In this paper, we concentrate on the thickness- 
driven SIT in this compound. The following sections will 
be organized as follows : first, section [III will detail the 
experimental procedures. Section IIIII will explain the 
finite-size scaling method we have used to analyze our 
results concerning the disorder-induced transition under 
non-zero magnetic fields and show that we have obtained 
surprisingly small critical exponents for the transition. 
Combining this analysis with our previously obtained re- 
sults [19], we infer the phase diagram for Nba^Sii-a^ (sec- 
tion [IVl). Finally, section IVl will provide a discussion on 
the interpretation of these sets of experiments on disor- 
dered superconducting thin films and on the domain of 
validity of the " dirty boson model" . 



II. EXPERIMENTAL PROCEDURE 

The NbSi films have been prepared under ultrahigh 
vacuum by e-beam co-deposition of Nb and Si. A series 
of four samples with stoichiometry Nbo.15Sio.85 and thick- 
nesses of 100, 50, 25 and 12.5 nm have been deposited 
onto sapphire substrates coated with a 50 nm-thick SiO 
under layer. The films were synthesized during a single 
run in order to have the samples' niobium concentra- 
tions as similar as possible. We also took special care 
over the control of the sample's parameters : the evap- 
oration was controlled in situ by a special set of piezo- 
electric quartz in order to precisely monitor the compo- 
sition and the thickness of the deposition. These two 
characteristics were then controlled ex situ by Ruther- 
ford Back Scattering (RBS) and the results well fitted 
with the in situ monitoring. Samples of the same sto- 
ichiometry with thicknesses down to 2.5 nm have been 
characterized by Atomic Force Microscopy and showed 
no sign of morphological granularity nor inhomogeneity. 
The superconducting transitions of these samples in zero 
magnetic field are a few tens of mK sharp and show no 
sign of reentrant behavior as usually observed for gran- 
ular systems. Besides, all samples showed the same re- 
sistivity at high temperature within 4%. All these argu- 
ments lead us to think that our samples are homogeneous 
in composition, non granular and only differ from one an- 
other by their thickness. This conclusion is corroborated 
by a TEM study [2l| showing that only Nba^Sii-a^ alloys 
annealed at 500C present Nb-rich clusters. The electri- 
cal characteristics of the four films were measured down 
to 150 mK using a dilution refrigerator. A perpendic- 
ular magnetic field could be applied and was made to 
vary between 5 and 11 kOe. Resistance measurements 
were performed using a standard AC lock-in detection 
technique operated at 23 Hz. A current of 100 nA was 
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TABLE I: Relevant parameters for our samples : the thickness 
d, the superconducting transition temperature Tco, the BCS 
coherence length ^o, the effective coherence length ^eff and 
the dephasing length L$ at 0.3 K. 



d [nm] 


Tco [mK] 


^0 [/im] 


^eff [nm] 


L$(0.3 K) [nm] 


12.5 


213 


12.8 


58.2 


50 


25 


347 


7.9 


45.7 


64 


50 


480 


5.7 


38.9 


75 


100 


530 


5.2 


37.1 


79 



applied to the sample, which is within the linear regime 
of the I-V characteristics for the considered films. All 
electrical leads were filtered from RF at room tempera- 
ture. 



III. (i-INDUCED TRANSITION 

Before describing the renormalization procedure we 
have used and the results thus obtained, let us establish 
the dimensionality of our samples. In our system, the 
mean free path / is of the order of the interactomic dis- 
tance : / :^ 2.65 A [22[ and hence much smaller than the 
superconducting coherence length given by the BCS 
theory (^o = 0.18^^ where vp is the Fermi velocity 
estimated to be 2x10^' cm.s"^ Q). In the "dirty" limit 
the effective coherence length of the system is given by 
^e// = We also have to consider the dephasing 

length which acts near the SIT as a cut-off length due to 
the finite temperature [H, 0, [ijj : = ^ where 
rrie is the mass of the electron. The smallest length be- 
tween L$ and ^eff hence determines the dimensionality 
of the film. The different relevant lengths are given in 
table H The films with thicknesses ranging from 12.5 to 
50 nm can be considered to be 2D, whereas the 100 nm 
film is 3D. In the renormalization procedure we shall fo- 
cus on the 2D films, so that the resistances mentioned 
below are sheet resistances. Let us also note that, in 
what follows, we used the usual convention found in the 
SIT-related literature [24]: the term "superconducting" 
applies to curves that have a positive Temperature Co- 
efficient of Resistance (TCR : ^), and, by contrast, we 
shall label as "insulating" all curves having a negative 
TCR. 

As shown by the rarity of experimental data concerning 
the thickness-induced SIT, it is difficult to obtain a series 
of samples that are identical except for their thickness : 
unlike the magnetic field, d cannot be tuned continuously. 
We have therefore developed an analysis method which 
enables us to interpolate the system's transport behavior 
between the discrete values of d we experimentally have 
access to. 

All four samples were superconducting at zero mag- 
netic field (insert of figure [1]) and were progressively 
tuned through the transition by a finite H. For each value 
of all four samples were studied (figure [1]) and the di- 
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FIG. 1: Resistance per square as function of temperature for 
if = 6.8 kOe. The curves for all four samples are represented. 
For this particular value of the magnetic field, the 25 nm, 50 
nm and 100 nm-thick films are superconducting, whereas the 
12.5 nm-thick film is insulating. Inset : The same data at 
zero magnetic field. 
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FIG. 2: Resistance per square as function of sample thickness 
for H = 6.8 kOe. The curves are represented for 16 different 
values of the temperature between 168 and 831 mK. Insert : 
the same data are shown around the crossing point at about 
dc = 23 nm for four particular temperatures : T = 168, 186, 
239 and 831 mK. This crossing point is interpreted as the 
signature of a QPT. 

agram (R^d) traced for different temperatures presents 
a crossing point (figure [2]). This is the signature of the 
QPT [10] and allows us to estimate the critical thick- 
ness dc associated to the magnetic field H. We repeat 
this process for all values of if, obtaining a collection of 
critical parameters couples {dc^H). 

When H is fixed, the thickness-induced transition is 
solely governed by the distance to the transition 5d = 
^^^^""^ . If these (i-driven transitions all belong to the same 
universality class, independent of the particular value of 
ii", the only relevant parameter for the scaling of all our 
data is the value of 5d = This means that all 

curves i?(^^^^^^,T) should collapse on two universal 
curves. Note that the renormalized quantity we consider 
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FIG. 3: Renormalization of the resistance R for the critical 
exponents — 0.4 for the 25 nm-thick sample. Each color is 
affected to a particular value of 5d- Left insert : determination 
of the critical exponents product by the derivation method. 
Right insert : determination of the critical exponents product 
by the t(T) minimization method. 



is R and not ^ as in 1(1 for we do not find an universal 
n — 

critical resistance For each individual sample, this 
means that by tuning i7, dc is made to vary and so does 
5d' In other words, the thickness d being fixed, the criti- 
cal thickness dc is changed via the magnetic field. Since 
the only relevant parameter for the scaling is the distance 
5d to the transition, this situation is ultimately equivalent 
to having a fixed critical thickness and variable sample 
thicknesses (as in [10] for example). 

For each sample, the results were analyzed using two 
independent scaling methods [7, 10]. First, for the deriva- 
tive method, we plot -^^I^^^q ^ RcT ^ f'{0) as func- 
tion of ^ which, in a log-log diagram, gives a straight 
line of slope — (left insert figure [3|). The second 
method consists in numerically finding t{T) such that 
R{6d,t{T)) = Rcf{Sdt{T)) and that t{T) yields the best 
collapse between the data measured at the temperature 
T and the data measured at our lowest temperature (150 

mK). To obey the scaling law (equation [1]) , t{T) should 
_ 1 

be of the form T and we can hence infer the value 
of VdZ (right insert figure [3]). 

For all 2D samples, we obtained a product of critical 
exponents of i/^z = 0.4 dz 0.15 . We can check this 
value of the exponents product by plotting R as function 

of * T ^ (figure [3|) for the 25 nm-thick sample. All 
data superimpose nicely in the ranges 0.16 < T < 0.35 
K and | Sd \< 1, forming two curves only, one repre- 
senting the superconducting behavior and the other the 
insulating side of the transition. \ Sd \ = I still exhibits 
a critical behavior since the corresponding data collapse 
on the same curves. It is quite surprising that the scaling 
continues to work that far from the critical point. The 
analysis performed on the 12.5 nm and the 50 nm-thick 
samples gave the same value of the product UdZ within 
the uncertainty. 
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FIG. 4: Renormalization of the renormalized resistance 
for the critical exponents UdZ = 0.4 for the 12.5 (triangles), 
25 (circles) and 50 (squares) nm-thick samples. 



This far, we have only considered the renormalization 
of the resistance for one sample at a time. In order to 
compare the critical behavior of the different samples, 
we have to take into account their different normal resis- 
tances. We therefore have to compare the quantity 
where Rn is the resistance taken at high temperature, 
typically at IK. This procedure is not usual in the lit- 
erature and directly derives from the fact that, in our 
experiment, Rc is not universal and varies over one order 
of magnitude (see section V). The scaling of then has 
no significance [3j. 

We then looked for a critical exponent product that 
allowed all curves from all samples to collapse. For each 
sample, we adjusted the non-universal parameter a of 
equation ([T]) for the curves to superimpose. We found 
Q^i2.5nm = 1-9, <^25nm = 0.9, and a^onm = 0.5 for a prod- 
uct of UdZ = 0.4 dz 0.1. The corresponding criteria for 
the renormalization are then very clearly defined : i. the 
magnetic field was made to vary between 5.1 and 10.5 
kOe by increments of 0.1 kOe ; all critical points (dc^Hc) 
corresponding to these fields have been taken into ac- 
count ; ii. the only constraint on the distance to the 
transition is S < 0.8 ; iii. 0.17 < T <0.39 K. The result 
of the renormalization is given figured! This graph is par- 
ticularly remarkable : even if our samples have normal 
resistances varying by nearly one order of magnitude, the 
corresponding resistances all collapse on a single renor- 
malization plot. 



IV. PHASE DIAGRAM 

The renormalization method has enabled us to mea- 
sure a number of critical parameters couples (Hc^dc) al- 
though we only had four different samples. We can hence 
draw part of the phase diagram for Nb^^Sii-a; thin films 
(figure [5]). The line formed by the critical points sepa- 
rates an insulating region at high fields and small thick- 
nesses form a superconducting region at low field and 
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FIG. 5: Phase diagram for a-NbisSiss in the (He, dc) plane. 
The open symbols were obtained from the thickness-tuned 
SIT, whereas the full symbols were obtained in [19] for the 
magnetic-field tuned transition. 



large thicknesses. Of course these critical points coincide 
with those determined from the magnetic field-induced 
SIT [19]. As for a-Bi [10|, depending on the parameter 
tuned to cross this line, the critical exponents product 
found is different : = 0.7 when the field is varied, 

whereas a variation of the sample's thickness gives = 
0.4. We thus confirm that these two SITs belong to two 
separate universality classes. 



V. DISCUSSION 

First let us comment on the value found for the critical 
exponents product. For a-NbSi in a thickness-induced 
SIT, we have found = 0.4. This value is surpris- 
ing when compared to other critical exponents found 
by other groups, thin a-Bi films for instance, for which 
UdZ = 1.4. At this point we do not have any clear expla- 
nation for this important difference. However = 1.4 
is close to what is predicted for classical 2D percola- 
tive systems {ud = 4/3) and a-Bi thin films present a 
thickness-induced SIT for very shallow thicknesses (a few 
angstroms, 20 A at most). In this sense also, our system 
is particularly interesting since it allows 2D samples to 
experience a thickness-driven SIT at reasonable thick- 
nesses where the roughness of the film, the surface state 
of the substrate or the microscopic details of the film's 
growth should not be important factors. UdZ = 0.4 is also 
surprisingly small considering the theoretical predictions 
that have been made to this day [1]. Although the exact 
value of this product might be affected by the uncertainty 
on the determination of the exponents (±0.1) and by the 
small number of samples we have, at any rate, we can 
confidently say that UdZ < 1 which is inconsistent with 
the "dirty boson model". If we assume that ^ = 1, the 
consequence of this is that Ud < Many authors have 
pointed to the fact that this violates the so-called "Harris 
criterion" {v > 2/d) [25], however this criterion is valid 



FIG. 6: Critical resistance as function of the critical field for 
a-NbSi films. 



for small disorder and since our system consists in amor- 
phous films in which the mean free path is of the order of 
the inter-atomic distance, it is not all that shocking that 
the value found for the localization length exponent does 
not obey this inequality [26]. 

Another point that has much been discussed related to 
the " dirty boson model" is the value of the critical sheet 
resistance. In this set of experiments, we show that Rc 
varies over a large range when either the magnetic field 
or the thickness is varied (figure [6|). 

Until now we have analyzed our results by comparing 
them to the "dirty boson model". Although the renor- 
malization procedure works remarkably well for all sys- 
tems studied to date - which means the SIT is indeed 
a QPT [1] -, two important predictions of this model 
{v > 2/d and Rc = h/Ae^) are not verified by a-NbSi 
thin films as well as in other systems (a-Bi [10], a-Be [9|, 
NdCeCuO [27.], MoGe [28], InOx [13, 16], TiN [12.], MoSi 
[8]...). One might therefore put this model int o q ues- 
tion. Tunneling effect experiments suggest ['2^, l3d Isij 
that, for homogeneous systems, amplitude fluctuations 
of the order parameter play a role even in the vicinity of 
the SIT : when the films' thickness decreases, the super- 
conducting gap A and the critical temperature decrease 
together, monotically, such that 2A/Tc ^ constant. In 
this picture, near the SIT, the amplitude of the supercon- 
ducting order parameter can become very small, whereas 
an essential point in the "dirty boson model" is that its 
amplitude is finite near the transition. The same stud- 
ies show that, even in the "superconducting" - in the 
previously-defined sense of the TCR - region, the one 
particle density of state is not zero, meaning that there 
are normal excitations coming from electrons that are 
not involved in any Cooper pair. This would mean that 
amplitude fluctuations of the system must be taken into 
account for a correct description of the transition, which 
is not the case in M.P.A. Fisher's model. The suggestion 
by some authors that other phase (s) may be involved in 
between the superconducting and the insulating regimes 
is particularly interesting. Some have suggested a vortex- 
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FIG. 7: Resistance as function of the temperature for the 50 
nm-thick sample dX H — 7.9 kOe. Over one decade variation 
in temperature, the film's resistance only varies within 3.5 Q 
(0.5% in relative value), which is our experimental uncertainty 
in this range. 



which vhz = 0.65. These two SITs therefore belong to 
two different universality classes. However the very small 
value of cannot be explained by the existing models 
for this transition. Further experimental investigations 
are needed to understand the growing discrepancies be- 
tween the various experimental results and between these 
results and the theory. 
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liquid phase j32|, [s^ which has recently [34] been linked to 
the problem of anomalous Nernst effect in the cuprates. 
As recent measurements on amorphous superconductors 
have shown [35, 36, 37], Nernst effect is a very sensitive 
probe of amplitude fluctuations |^ and phase fluc- 
tuations [37^ of superconducting order parameter. These 
last works suggest that measurements of the Nernst ef- 
fect should be a relevant probe to test the existence of 
this vortex liquid phase. However, there has not been 
clear predictions on how the thickness variation should 
affect this phase. Also very appealing is the suggestion 
that there is a bosonic metallic phase, such as the Bose 
metal [l8|, involved. This hypothesis is very interesting, 
in particular when one takes a close look at the resis- 
tive behavior of our films. Indeed, at low temperatures, 
the resistance of some samples seem to saturate at a fi- 
nite value (figure [7j), displaying a large temperature range 
where the resistance is independent of the temperature. 
However a study at lower temperatures should be under- 
taken to confirm this tendency. Let us restate that the 
qualification of insulating or superconducting have been 
arbitrarily attributed to < (resp. > 0) curves 
without any other ground than the assumption made by 
the "dirty boson model" that only these two phases ex- 
isted. All these arguments (the amplitude fluctuations 
of the order parameter, a possible fermionic channel, the 
suggestion of a Bose metal...) plead in favor of a recon- 
sideration of the " dirty boson model" and further exper- 
imental investigations of these systems. 

In conclusion, we have studied the thickness-induced 
SIT in the presence of a perpendicular magnetic field on 
a-NbisSigs thin films of thicknesses ranging from 12.5 to 
100 nm. We have found the signature of a QPT when 
the sample thickness is lowered. The corresponding crit- 
ical exponents product is — 0-4 ±0.1. This value is 
different from the one found in the analysis of the mag- 
netic field-induced transition in the same compound for 
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